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1.  INTRODUCTION 


When  designing  adhesively  bonded  fibre  composite  repairs  for  metallic  or  com¬ 
posite  structures,  ttvo  of  the  main  design  requirements  are: 

(a)  the  adhesive  bond  must  not  fail 

(b)  the  composite  repair,  often  referred  to  as  a  patch  [1],  must  not  delaminate. 
These  failure  modes  involve  matrix  dominated  failures  and  require  a  consistent  and 
>alid  failure  criterion  for  use  in  the  design  phase  of  any  repair. 

A  similar  situation  occurs  in  the  compressive  failure  of  impact  damaged  lami¬ 
nates.  This  failure  process  is  often  matrix  dominated  [2]  and  has  received  consid¬ 
erable  attention  in  recent  years.  A  review  of  the  current  status  of  bonded  repairs 
and  the  testing  and  analysis  methodologies  for  impact  damaged  laminates  is  given 
in  [3,4]  respectively. 

The  pvrpose  of  this  present  paper  is  to  outline  the  main  approaches  currently 
used  for  the  analysis  of  matrix  dominated,  failures  and,  where  possible,  to  explain 
the  relationship  between  them. 

2.  FAILURE  THEORIES 

There  are  a  large  number  of  failure  theories  presently  used  to  characterize  the 
matrix  controlled  failure  of  composites  and  adhesives.  The  theories  that  will  be 
discussed  in  the  present  paper  are, 

(1)  the  tensor  polynomial  failure  criteria,  i.e.  Tsai  -  Hill  and  Tsai  -  Wu, 

(2)  the  energy  density  failure  hypothesis,  and 

(3)  the  energy  release  rate  approach. 

The  first  two  approaches  differ  from  the  latter  in  that  they  involve  the  stresses, 
and/or  strains,  evaluated  at  a  characteristic  distance  rc  in  front  of  the  stress  concen¬ 
trator.  The  energy  release  rate  approach  does  not  involve  a  length  scale,  which  is 
a  material  parameter,  but  instead  requires  a  series  of  additional  tests  to  determine 
the  mode  I,  II  and  III  energy  release  rates. 

2.1  The  Tensor  Polynomial  Failure  Criteria 

In  recent  years  numerous  forms  of  the  polynomial  failure  criteria  have  been 
proposed.  The  most  generally  accepted  of  these  can  be  expressed  in  the  form 


F;<r,  +  Fij<?i<Tj  +  Fijk<Ji<TjOk  =1  ( 1) 

where  F,,  Fij  and  Fy*  are  the  constants  which  must  be  experimentally  determined, 
see  [5]. 

When  a  body  contains  a  stress  concentrator,  such  as  a  matrix  crack,  then  Wu  [G] 
stated  that  this  failure  criterion  should  not  be  evaluated  at  the  stress  concentrator 
but  at  a  distance  »y  from  it.  Indeed  it  was  shown  that  ly  is  a  material  constant 
which  can  be  obtained  experimentally. 

In  the  case  of  the  matrix  dominated  failure  of  a  fibre  composite  laminate  Hahn 
and  Tsai  [7]  subsequently  proposed  that  the  failure  criterion 

Fi*i  +  FijViOj  =1  («  =  2, ...,  6)  (2) 
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replaces  equation  (1).  The  main  difference  here  lies  in  the  neglect  of  the  contribution 
of  the  fibre  stresses,  i.e.  9\.  Equation  (2)  is  still  evaluated  ».t  a  distance  re  from  a 
stress  concentrator. 

This  theory  has  been  successfully  used  by  many  reseachers  and  has  been  applied 
to  the  problem  of  delamination  [8]  as  well  as  to  the  fracture  of  centre* notched  panels 
[5]. 

If  the  problem  under  consideration  does  not  involve  shear  stresses  then  equation 
(1)  reduces  to 

fjffi  +  =  1  (3) 

which  is  similar  to  the  point  stress  failure  criterion  [9]  and  which  has  been  extensively 
used  to  characterise  matrix  dominated  failure  [10,11,12].  In  references  using  the 
point  stress  formulation  the  term  rc  is  often  denoted  as  d0. 

For  a  general  problem  it  is  necessary  to  determine  the  failure  location  as  well  as 
the  failure  load.  One  method  for  locating  the  point  of  failure  initiation  was  proposed 
by  Wu  [5].  This  approach  considers  both  the  strength  vector  F  to  the  failure  surface, 
as  given  by  equation  (1),  and  the  stress  vector  S  where, 

\s\  *  (5) 

Failure  is  then  said  to  occur  at  the  location  where 


5|r.r,  =  F\r..' 


Details  of  this  hypothesis  are  given  in  [5]. 

When  the  material  is  isotropic  the  strength  function  F  reduces  to  the  modified 
von  Mises  yield  criterion 


W,.  =  AdV  +  B  (7) 

where  IV'j  is  the  distortional  energy,  A  and  B  are  constants  and  dV  is  the  change  in 
volume  per  unit  volume  i.e., 


dV  =  (<Ti  +  <rj  +  <tj)/3/\  (8) 

where  K  is  the  bulk  modulus.  This  form  of  the  yield  criterion  has  been  widely  used 
for  the  analysis  of  polymers  and  adhesives  [13, 1*1]. 

2.2  Energy  Density  Theory 

The  concept  of  a  critical  energy  density  has  been  extensively  used  in  the  design 
and  failure  analysis  of  composites  and  adhesively  bomb'd  repairs.  Applications  of 
this  theory  range  from 

(a)  fatigue  life  evaluation  of  damaged  laminates  [15] 

f  b)  evaluation  of  the  effects  of  combined  mechanical  and  enviromv.ital  loads  on 
hole  elongation  and  subsequent  failure  [16] 

(c)  estimating  the  residual  strength  of  impact  damaged  laminates  [17,18]  and 
laminates  containing  edge  delaminations  [19] 

(d)  characterizing  the  fracture  behaviour  of  fibre  composite  laminates  [20.21,22]. 
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It  Km  alto  been  used  to  design  adhesively  bonded  repair*  for  cracked  metallic  com¬ 
ponents  [23]  and  was  used  to  design  a  boron/epoxy  reinforcement  to  the  wing  pivot 
fitting  ot  F-l  11  aircraft  in  service  with  the  Royal  Australian  AirForce  [24] 

This  theory  also  involves  a  length  scale  re  and  states  that,  for  matrix  dominated 
failures,  failure  initiates  in  the  direction  6o  which,  at  r  *  re,  coincides  with  tha 
local  minimum  of  the  energy  in  the  matrix,  defined  as  Wt.  Failure  occurs  when  this 
value  reaches  a  critical  value  We.  Mathematically  this  can  be  written  as 


where  Wt  is  given  by 


0) 

W,  =t  We  at  ©  =  r  =  rc 

(10) 

VV«  =  l/2ffy<y  -  Wj 

(ID 

where  Wj  is  the  energy  in  the  fibres.  When  applying  this  hypothesis  to  the  failure 
of  adhesively  bonded  joints  Wj  is  considered  to  be  sero. 

As  in  the  tensor  polynomial  failure  criterion  this  approach  involves 

(1)  the  use  of  a  characteristic  dimension  re 

(2)  the  neglect  of  the  contribution  due  to  the  fibres. 

Now  consider  the  failure  of  a  centre  notched  panel  where  the  fibres  are  parallel 
to  the  direction  of  the  crack.  It  follows  from  the  solution  to  this  problem  [2P1  that 
the  value  of  re  used  :n  the  energy  density  theory  is  related  to  that  used  in  the  tensor 
failure  theory,  which  we  will  denote  as  r* ,  by  the  expression 


*c  =  r‘U  -  (12) 

This  allows  the  test  methodologies  developed  for  the  measurement  of  r‘  to  be  used 
to  evaluate  rc.  As  a  first  estimate  it  is  often  convenient  to  use  rc  =•  r*. 

The  critical  energy  VV’f  is  a  function  of  d  V  and  the  level  of  local  constraint,  where 


J,...,  1  _  "lj,,,,  1  "13  "33.  .  1  "13 

^  '  p  p  r  ^  ^1'  r  p  r  ^  r  p 

£|1  ^11  *11  C;j  fin  fijl  fijj  fijj  Cl, 

This  dependence  is  shown  in  Table  2  where  the  strengths  and  moduli  used  are  taken 
from  Tsai  [26]  and  are  given  in  Table  1. 


"  ible  It  COMPOSITE  PROPERTIES 


MATERIAL 

"u 

l 

1 

T300/5208 

181 

10.3 

0.28 

7.17 

1.44 

2 

B4/5502 

>04 

18.5 

0.23 

5.59 

3.31 

3 

AS/3501 

138 

8.96 

0.3 

7.1 

1.26 

Note:  values  in  Table  are  in  GPa 
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Tabic  2:  DESIGN  ALLOWABLES  -  ORTHOTROPIC  APPROACH 


TENSION  COMPRESSION  SHEAR 


MATERIAL 

w 

d? 

w 

d? 

*.» 

w 

d.~ 

MPa 

MPa' 

MPa 

MPa 

MPa' 

MPa 

MPa 

MPa' 

MPa 

1 

40 

800 

27.4 

246 

30258 

-168 

68 

3329 

0 

2 

61 

1860 

41.4 

202 

20402 

-137.2 

67 

7429 

0 

3 

51.7 

1336 

35.2 

206 

21218 

-140.2 

93 

5449 

0 

where:  W  m  W  x  Eu  dV  =  dV  x  E„ 
Note:  is  taken  as  0.3 


The  failure  envelope  predicted  for  AS/3501,  using  both  the  energy  density  and 
Tsai  -  Hill  theories,  is  shown  in  Figure  1.  Both  theories  predict  essentially  the 
same  behaviour.  Figure  2  shows  the  critical  available  energy  density  for  the  three 
composite  materials  considered  in  Table  1 . 

As  mentioned  in  section  2.1,  the  Tsai  -  Hill  theory  is  based  on  a  particular  gen¬ 
eralisation  of  the  modified  von  Mises  criterion  for  isotropic  material.  An  alternative 
approach,  for  isotropic  materials,  is  to  add  the  energy  due  to  a  change  in  volume, 
denoted  by  to  both  sides  of  equation  (7)  and  then  noting  that 


EdV 1 

w  = 

"  6(1-2*) 


(H) 


and  that  the  total  energy  density  VV  =  W'^  +  IV,,  we  obtain  an  alternative  form  of 
the  modified  von  Mises  yield  criterion,  vii: 


W  =  Ad\f  +  B  +  EdV'/6(l  -  2*)  =  F{dY)  (15) 

For  any  given  material  the  function  F(dV)  should  be  experimentally  determined. 
This  equation  has  been  used  in  [24]  to  design  bonded  repairs  to  damaged  metallic 
components. 

If  we  replace  F(dV)  by  VVf  we  see  that  both  equations  (1)  and  (15)  represent 
possible  alternative  extensions  of 'he  modified  von  Mises  yield  criterion. 

The  concept  of  a  critical  energy  density  forms  the  l>asis  for  another  approach  to 
the  design  of  adhesively  bonded  joints,  see  [13],  In  this  approach  the  true  stress  - 
strain  curve  is  replaced  by  a  simpler  curve  with  a  constant  pt*t  yield  slope.  This 
idealized  curve  is  chosen  to  have  the  same  energy  density  as  the  true  stress  -  strain 
curve.  It  is  also  assumed  that,  for  asymmetric  joint,  the  only  stresses  in  the  adhesive 
are  shear  stresses.  These  approximations  allow  a  closed  form  analytical  solution  to 
be  obtained.  In  this  approach  the  shear  strain  in  the  adhesive  is  chosen  as  the  critical 
design  variable.  However,  because  of  the  choice  of  the  idealized  stress  -  strain  curve 
there  is  a  one  to  ore  relationship  between  shear  strain  and  energy  density. 
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1.3  Etscrgy  R*h—  Rat*  AppntdtM 

For  mode  I  self  similar  crack  growth  of  a  through  crack  in  an  iaotropie  body  in 
the  absence  of  body  forces  Atiuri  [27]  shovred  that  the  energy  release  rate  G  can  be 
written  as, 


G  m  Um  /  [Win  -  (16) 

r-OJr. 

where  n*  is  the  component  of  the  unit  normal  to  the  path  in  the  z,  direction,  tj  is  the 
traction  tensor  defined  as  f j  ■  W  is  the  strain  energy,  «,  is  the  dispiace.nent 
tensor  and  f,  is  a  vanishing  small  path  surrounding  the  crack  tip. 

It  has  also  been  shown  [27]  that  for  a  general  problem  'his  integral  will  not  equal 
the  value  of  G  obtained  from  the  movement  of  load  points  method,  i.e., 


r.p 

°  *  2B  d.x 


(17) 


where  P  is  the  applied  load,  C  «  S/P,  6  is  the  movement  of  the  load  points,  a  is 
the  crack  length  and  B  is  the  thickness  of  the  structure. 

This  formulation  was  extended  to  composite  materials  by  Jones  [28).  In  this 
case  G  is  given  by  the  same  expression  as  in  (16)  with  W  replaced  by  the  energy 
density  for  a  composite  material,  i.e., 


W  =  -  Ha'ijiT  -  To)  -  4w(M  -  Mo)  +  C,(7\  M)  (18) 

where  M  :s  the  mass  flux  of  moisture  per  unit  volume,  fa  and  ^  are  the  coefficients 
related  to  the  thermal  and  moisture  expansion  coefficients  of  the  body  respectively. 
To  and  Mo  are  the  reference  tempciature  and  moisture  content  respectively,  Cytj  is 
the  stifTness  tensor  and  C\  is  *  function  of  the  temperature  and  moisture. 

This  formulation  assumes  that  the  terms  in  the  integrand  which  are  functions 
of  temperature  and  moL.ure  alone  integrate  to  sero.  As  shown  in  [29]  the  existence 
of  body  forces,  which  for  matrix  dominated  problems  arise  due  to  closure  [30],  again 
result  in  compliance  measurements  giving  inappropiate  measurements  of  G- 

In  three  dimensions  the  integral  on  the  right-hand  side  of  equation  (16)  no 
longer  equals  G  and  is  referred  to  as  T*,  see  [27].  Approaches  based  on  energy 
release  rate  considerations  have  been  widely  used  in  characterising  matrix  failure 
[31,32,8,30].  However,  as  mentioned  above,  it  must  be  emphasised  that  compliuitce 
measurements  do  not  always  characterise  near-tip  behaviour.  This  is  particularly 
true  in  the  case  of  combined  mechanical  and  environmental  loading,  see  [28,33].  In 
this  case  the  equation  governing  the  coupling  between  mechanical  deformation  and 
thermal  energy  is  given  by, 


-hi,i  +  pfM  =  pCvT  -  (<v>-^  -  ~  7b')  -  “*y(M  -  M0)))Ti*(19) 

where  hi  is  the  heat  flux  tensor,  q  is  the  heat  absorbed  when  1  gram  of  moisture 
is  absorbed  by  1  gram  of  material  and  p  is  the  density,  see  [28].  Experimental 
confirmation  of  this  formulation  is  given  by  Wong  [34]  for  the  special  case  when  the 
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nwtttiil  is  isotropic.  If  the  mater  al  is  cyclic  ally  stressed  at  a  frequency  such  that 
adiabatic  conditions  hold  then 


PfSl -pi\T  -  -  To))  -  -  U»)))T:ij  (20) 


i) 


Let  u$  now  consider  a  bar  subjected  to  uniaxial  strain. 


cjj  m  At  strict  c  ■  0  if  »',>  ft  2  (21) 

substituting  into  (20)  gives, 

pqil  -  pCvT  =  (d}i*A<ro*wf  +  r  sin2*i1)T  (22) 

*  2  oT 

and  integrating  with  respect  to  t  we  obtain: 

"A M  -  pCvAT  =  ( A, Acs.W  -  |  A<*(  1  -  cos2ud))7o  (23) 

4  a/ 

This  shows  that  although  the  load  is  applied  at  a  single  frequency  w  the  temper¬ 
ature  and  moisture  fields  have  a  response  at  both  w  and  2w  and  that  the  amplitude 
of  the  response  at  frequency  2u<  is  proportional  to  the  square  of  the  local  strain. 
Hence  for  a  problem  involving  enviromental  effects  the  load  point  behaviour,  which 
only  reflects  the  behaviour  at  the  applied  frequency  w,  cannot  provide  information 
on  the  component  of  the  near-tip  field  which  is  responding  at  a  frequency  of  2<*>. 
Furthermore  the  near-tip  component  has  an  amplitude  proportional  to  At1  and  can 
thus  be  expected  to  be  very  large.  Near  the  tip  the  large  strain  field,  theoretically 
infinite,  means  that  adiabatic  conditions  are  never  realised  so  that  the  fatigue  be¬ 
haviour  will  depend  on  test  frequency.  Indeed  this  has  been  confirmed  by  a  series 
of  recent  laboratory  tests  [35]. 

The  existence  ot  a  response  in  the  temperature  field  at  frequencies  w  and  2w  has 
been  confirmed  experimentally  in  [36]  using  thermal  emission  techniques. 

U  must  be  stressed  that  energy  release  rate  approaches  are  best  used  when 
growth  is  co- planar  and  self-similar.  In  other  circumstances  they  must  be  used  with 
caution  [28,37].  An  alternative  approach,  which  is  related  to  the  energy  release  rate 
method,  was  suggested  by  Watanabe  [38],  This  method  uses  the  first  term  in  the 
expression  for  7V,  which  we  will  define  as  V/,,  to  characterise  the  damage,  vis  : 

W,  m  iim  i  Wdy  (24) 

*— °/r. 

This  formulation  has  not,  as  yet,  been  widely  used  but  when  applied  to  the  failure 
of  a  damaged  fastener  hole  under  both  thermal  and  mechanical  loads  appeared  to 
correlate  with  the  observed  failure  load,  see  [28],  The  parameter  W(  coincides  with 
T*  if  the  delamination  has  a  blunt  tip.  Further  work  is  required  oefore  this  can  be 
considered  a  valid  fracture  parameter. 

A  detailed  review  of  current  energy  release  rate  test  methodologies  is  given  in 
[2],  whilst  certain  problem  areas  arising  when  using  energy  release  rate  approaches 
are  highlighted  in  [37]. 
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Il  i»  interesting  io  note  that  the  paint  stress  and  average  stir**  criteria  run  also 
be  used  la  obtain  estimates  of  the  energy  release  rale  using  either  doitlJe  cantilever 
or  centre-notched  tent  specimens,  nee  [10,1 1).  The  values  «>f  C*V  obtained  in  this  way 
tend  to  be  higher  than  thoae  obtained  front  compliance  measurement*.  see  [2, 10] 

It  ia  well  known  that  for  metals  the  critical  sirens  intensity  factor  AV  is  strongly 
dependent  on  the  level  of  local  restraint.  Indeed  recent  tests  at  the  t’niversity  of 
Melbourne  [39]  have  shown  that  A't  for  an  6inm  thick  2024  TO  altintiniunt  alloy  is  ap¬ 
proximately  52  MPa^m  whilst  its  plane  strai-  alue  is  approximately  30  K!  Pr  >/m. 
This  corresponds  to  an  eight  fold  increase  in  f»,.  However,  this  effect  is  often  over¬ 
looked  for  composites.  Similarly  io  a  structure  the  level  of  rest  raint  is  variable  ami 
that  experienced  by  an  edge  delamination  may  be  quite  different  to  that  experienced 
by  a  delamination  surrounding  a  fastener  or  a  joint.  It  is  clear  that  a  methodology 
is  required  to  quantify  this  level  of  restraint.  One  such  methodology  based  on  the 
local  change  of  surface  and  volume  energies  is  presented  in  [40]. 

In  the  case  of  a  crack  lying  entirely  within  ui  adhesive  layer  the  energy  release 
rate  G  can  be  related  to  the  energy  density  W  by, 


G 


2(1  -y) 
(1-2*) 


r<w 


(25) 


However,  as  seen  in  section  2.2,  the  critical  value  cf  W  may  be  dependent  on 
dV.  This  infers  that  Gc  may  be  a  function  of  dV.  In  the  case  of  a  delamination  in  a 
composite  it  is  well  known  that  Gje  «  much  larger  than  G\t,  see  [2].  It  is  therefore 
tempting  to  conjecture  that, 


„  „  dV  ,,  <fl\  _ 

c*  u5vii,..  +(  5vT)u"  ** 


(26) 


where  dV \  is  the  value  of  dV,  at  r  =  rc,  which  occurs  when  Ge  --  G le.  However,  this 
expression  lacks  a  physical  basis  and  its  usefulness  in  describing  mode  III  failure 
is  questionable.  Tb  date  the  majority  of  methods  used  to  determine  Ge  for  mixed 
mode  failure  an  empirical  and  require  a  knowledge  of  (he  individual  mode  I,  II  and 
III  energy  release  rates  as  well  as  their  respective  failure  energies. 


3.  MATERIAL  NONLINEARITIES 


In  the  previous  sections  it  has  been  assumed  that,  up  to  failure,  the  material  is 
behaving  elastically.  However  a  recent  investigation  indicated  that,  prior  to  failure, 
there  was  significant  dissipated  energy  [41].  To  confirm  this  a  series  of  tests  was 
performed  on  a  56-ply  XAS-914C  laminate  with  a  ply  configuration  of  [±45/0j]r,. 
Six  specimens  with  dimensions  approximately  315  by  6.0  by  4.8  mm  were  prepared 
and  loaded  in  an  500  kN  Inst,  on  testing  machine  at  a  loading  rate  of  600  kN/sec. 
This  rate  of  loading  was  chosen  in  an  attempt  to  achieve  adiabatic  conditions. 

For  each  specimen  the  temperature  of  a  centrally  located  point  on  the  sur¬ 
face  of  the  specimen  was  measured  using  an  infrared  detector.  A  typical  non- 
dimensionalixed  temperature  response  is  shown  in  Figute  3.  In  each  case  failure 
of  the  specimen  did  not  occur  near  the  point  at  which  the  temperature  was  mea¬ 
sured.  Details  ci  the  test  methodology  are  given  in  [42]. 

From  Figure  3  it  is  seen  that  the  measured  temperature  response  is  simila.'  to 
that  for  metals  [43].  There  is  an  initial  cooling  during  the  linear  regime,  which  for 
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imtik  is  described  by  Kelvin’s  law,  fallowed  by  extensive  heating.  The  amount  of 
heating  ie  a  direct  measure  of  the  dissipated  energy,  and  as  such  supports  the  results 
given  in  [41]. 

As  »  result  of  this  investigation  it  is  clear  that  farther  work  is  necessary  i mirier 
to  understand  the  role  of  dissipated  energy  in  the  failure  process. 

4.  CONCLUSION 

This  paper  has  outlined  several  methods  currently  used  for  analysis  of  the  matrix 
dominated  failure  in  composite  materials  and  adhesive!)-  bonded  joints.  Particular 
attention  has  been  given  to  the  relationship  between  these  approaches. 

Attention  has  also  been  focused  on  the  role  of  dissipated  energy  in  the  failure 
process.  This  is  a  topic  which  is  not,  as  yet,  fully  understood  and  with  the  advent 
of  more  ductile  matrix  materials  requires  greater  attention. 
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